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A calculation of site-bond percolation thresholds in many lattices in two to five dimen- 
sions is presented. The line of threshold values has been parametrized in the literature, 
but we show here that there are strong deviations from the known approximate equa- 
^ ' tions. We propose an alternative parametrization that lies much closer to the numerical 

OO ' values. 
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^^ ' 1. Introduction 

^' 

On . Site-bond percolation is a natural generalization of pure site percolation and pure 

bond percolation. The generalization allows both sites and bonds to be randomly 
occupied, with probabilities ps and pb respectively, in the case of random percola- 
tion. In particular, ii Pb — Ps, site-bond percolation reduces to a pure site or pure 

'^ ' bond percolation on a decorated lattice. It is also possible to introduce a correlation 

G . between the occupation probabilities of adjacent sites and bonds, which is usually 

referred to as correlated percolation. Site-bond percolation has many applications 
in different fields. Foii^example, it was used to describe the sol-to-gel transition 
(gelation) of polymers,tl the spreading of an infection,H etc. 

^\f ' In site-bond percolation one looks for clusters of occupied sites, connected to 

S . each other by occupied bonds. When there is a cluster that connects one side of the 

system with the other, the system is said to percolate. First mentioned by Frisch 
and Hammersleyjlthe model was investigated at first by Hoshen by means of Monte 
Carlo simulations.Q Agrawal et al.l3 showed, using a series method, that the critical 
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exponents of nure site percolation are also valid for site-bond percolation. Nakanishi 
and ReynoldsQ confirmed this conclusion through a position-space renormalization 
group study. They also found that the boundary between percolation and non- 
percolation in the {pb,ps) plane is governed by a single fixed point. Hanimersleya 
proved for a partially directed graph a theorem. The theorem yields the inequality 

P{PsPb, 1) < P{Ps,Pb) < P{l,PsPb), (1) 

where P{ps,Pb) is the percolation probability (i.e. the probability that a single 
source site is connected to an infinite set of other sites). The inequality Eq. (|l|) 
gives reasonably sharp bounds for the mixed percolation probability. For a Bethe 
lattice (pcb = Pes = Pc, where Pcs is the threshold for pure site percolation, and pcb 
is the one for pure bond percolation) we can get 

PsPb = Pc (2) 

The shape of the boundary between percolation and non-percolation was subse- 
quently studied by Yanuka and Englman.13 They proposed the following equation 
for the 'critical curve' in the (j>b,Ps) plane. 

\0gps , }0gpb _ . fo\ 

log Pes log Pcb 

The equation satisfies the inequality Eq. (|l|) and gives the correct limit Eq. (0) for 
a Bethe lattice. The equation defines a line in the (pb^ps) plane, starting at {pcb, 1) 
and ending at (l,pcs), as is shown in Figure ll|. Yanuka and Englman showed plots 
of Eq. (y) together with simulation results for {pb,Ps) for the square, triangular, 
simple cubic and face centered cubic (fee) lattices. All their simulation results lie 
remarkably close to the line defined by Eq. (j^). It is not clear from their paper 
whether the data of e.g. the fee lattice, which lie visibly below the line given by 
Eq. (^, are accurate enough to conclude that there is a discrepancy. 

Tretyakov and InuiB studied the same critical curve, but now for directed perco- 
lation. Their simulation results were accurate enough to conclude that deviations 
in the range lO^'' — 10^^ occur. This raises the question whether such deviations 
also occur for 'regular' site-bond percolation. 

Therefore we have undertaken a systematic study of the threshold values (j)b,Ps) 
for random site-bond percolation in many lattices. The results, which are presented 
in Section , show many deviations from the Yanuka-Englman line of Eq. (^. The 
deviations are largest for lattices with a low ratio of a = Pcb/Pcs, i-e. lattices 
for which the bond percolation thresholds is much lower than the site percolation 
threshold. In Section we discuss other approximative equations for the critical 
line {pb,Ps)- 



2. Numerical results 



All numerical results were obtained with the programs described in Ref.tr For 
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Fig. 1. The percolation threshold pairs {ps,Pb) for many different lattices in various dimensions. 
The points are from the simulations (see Tables Ill4l|) , the lines are from Eq. (tsl) . 
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Table 1. The results for the site-bond percolation thresholds for several lattices in two dimensions. 
Error estimates concerning the last digit(s) are indicated between brackets. 





square 


trian 


gular 


Kagome 


p 


Pb{Ps=p) 


PsiPb^p) 


Pb{Ps=p) 


Ps{pb=p) 


PbiPs^p) 


PsiPb^p) 


0.95 


0.5334(3) 


0.6149(3) 


0.3732(3) 


0.5152(4) 


0.5615(4) 


0.6711(4) 


0.90 


0.5716(3) 


0.6396(4) 


0.4023(3) 


0.5331(4) 


0.6046(5) 


0.6914(5) 


0.85 


0.6154(3) 


0.6673(3) 


0.4368(3) 


0.5518(4) 


0.6556(6) 


0.7162(5) 


0.80 


0.6663(3) 


0.6980(3) 


0.4764(4) 


0.5736(4) 


0.7152(7) 


0.7428(5) 


0.75 


0.7257(4) 


0.7322(4) 


0.5251(4) 


0.5969(4) 


0.7894(9) 


0.7757(8) 


0.70 


0.7961(4) 


0.7709(3) 


0.5830(6) 


0.6231(3) 






0.65 


0.8807(4) 


0.8157(3) 


0.6561(6) 


0.6541(4) 






0.60 


0.9840(4) 


0.8671(3) 












octag 


;onal 


honeycomb 


star-dual 


P 


PbiPs^p) 


PsiPb^p) 


Pbij>s=p) 


Ps{Pb=p) 


PbiPs^p) 


PsiPb^p) 


0.95 


0.3472(4) 


0.5133(7) 


0.6926(5) 


0.7275(5) 


0.2796(4) 


0.5036(5) 


0.90 


0.3744(5) 


0.5276(9) 


0.7377(5) 


0.7610(5) 


0.3022(5) 


0.5101(5) 


0.85 


0.4087(5) 


0.5442(7) 


0.7890(5) 


0.7986(5) 


0.3296(4) 


0.5188(5) 


0.80 


0.4449(6) 


0.5620(7) 


0.8481(5) 


0.8401(5) 


0.3615(5) 


0.5285(5) 


0.75 


0.4912(7) 


0.5816(7) 






0.4006(4) 


0.5417(5) 


0.70 


0.5481(10) 


0.6052(7) 






0.4499(5) 


0.5572(5) 


0.65 


0.6213(24) 








0.5133(4) 


0.5763(5) 


0.60 










0.5968(5) 


0.5984(5) 




puzzle 


direct 


low 


-pb 


P 


PhiPs^p) 


Ps(Pb^p) 


Pb{ps=p) 


Ps{Pb=p) 


PbiPs^p) 


Ps(Pb=P) 


0.95 


0.5693(5) 


0.4495(4) 


0.5677(4) 


0.4433(4) 


0.2166(5) 


0.5003(5) 


0.90 


0.5901(5) 


0.4849(5) 


0.5877(4) 


0.4780(4) 






0.85 


0.6124(4) 


0.5247(5) 


0.6091(4) 


0.5183(5) 


0.2554(5) 


0.5048(5) 


0.80 


0.6371(5) 


0.5720(5) 


0.6343(4) 


0.5648(4) 






0.75 


0.6655(5) 


0.6283(5) 


0.6602(4) 


0.6199(5) 


0.3110(5) 


0.5156(5) 


0.70 


0.6976(4) 


0.6952(5) 


0.6913(4) 


0.6876(5) 






0.65 


0.7328(5) 


0.7766(5) 


0.7263(4) 


0.7692(5) 


0.4020(5) 


0.5334(5) 


0.60 






0.7674(5) 


0.8712(7) 






0.55 










0.5644(5) 




0.45 










0.6149(5) 




0.35 










0.6997(6) 





An Investigation of site-bond percolation 5 



Table 2. The results for the site-bond percolation thresholds for several lattices in three dimensions. 
Error estimates concerning the last digit(s) are indicated between brackets. 



P PbiPs^p) PsiPb^p) Pb{Ps 



fee 

-P) Ps{Pb-' 



--P) Pb{Ps 



hep 

--P) Ps 



(Pb^p) 



0.95 
0.90 
0.85 
0.80 
0.75 
0.70 
0.65 
0.60 
0.55 
0.50 
0.45 
0.40 
0.35 



0.2637(3) 
0.2806(3) 
0.2991(3) 
0.3208(3) 
0.3452(3) 
0.3735(3) 
0.4072(3) 
0.4475(3) 



0.3244(3) 
0.3385(3) 
0.3539(3) 
0.3716(3) 
0.3912(3) 
0.4139(3) 
0.4395(3) 
0.4694(3) 



0.1273(3) 
0.1355(3) 
0.1453(3) 
0.1560(3) 
0.1679(3) 
0.1823(3) 
0.1995(4) 
0.2201(4) 
0.2443(5) 
0.2769(5) 
0.3161(6) 
0.3676(8) 
0.4404(11) 



0.2049(5 
0.2124(4 
0.2202(4 
0.2291(5 
0.2391(5 
0.2515(5 
0.2644(5 
0.2796(6 
0.2974(6 
0.3190(6 
0.3451(6 
0.3759(6 
0.4171(7 



bee 



p Pb{ps=p) Ps{j)b=p) 



diamond 
Pb{Ps=p) Ps{pb=P 



0.1272(3 
0.1361(3 
0.1456(3 
0.1559(3 
0.1688(3 
0.1832(3 
0.2001(3 
0.2209(3 
0.2455(3 
0.2764(3 
0.3164(3 
0.3698(3 



2052(3 
2125(3 
2214(3 
2306(3 
2404(3 
2519(3 
2653(3 
2806(3 
2985(3 
3195(3 
3454(3 
3764(3 



0.95 
0.90 
0.85 
0.80 
0.75 
0.70 
0.65 
0.60 
0.55 
0.50 
0.45 
0.40 
0.35 



0.1909(3 
0.2030(3 
0.2164(4 
0.2324(4 
0.2498(4 
0.2711(4 
0.2953(5 
0.3249(6 
0.3603(6 
0.4036(7 
0.4578(9 
0.5300(10 
0.6281(13 



0.2545(5 
0.2656(5 
0.2777(5 
0.2900(5 
0.3043(5 
0.3210(6 
0.3408(6 
0.3633(6 
0.3896(7 
0.4199(7 
0.4562(7 
0.5043(7 
0.5639(7 



0.4118(4) 
0.4366(4) 
0.4656(4) 
0.4973(4) 
0.5345(5) 
0.5758(6) 
0.6253(6) 
0.6845(7) 
0.7555(8) 



0.4489(5 
0.4723(5 
0.4959(5 
0.5229(5 
0.5544(6 
0.5885(5 
0.6283(5 
0.6755(6 
0.7308(6 
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Table 3. The results for the site-bond percolation thresholds for the 4-diniensional simple cubic 
(sc) and fee lattices. Error estimates concerning the last digit(s) are indicated between brackets. 



Pb 



4d- 

(Ps = P) 



SC 
PsiPb =P) 



Pb{Ps 



4d- 

p) 



fee 

Ps{Pb 



P) 



0.95 


0.1701(6) 


0.2065(6) 


0.0520(03) 


0.0853(08) 


0.90 






0.0554(04) 


0.0892(09) 


0.85 


0.1900(6) 


0.2254(6) 


0.0584(04) 


0.0920(10) 


0.80 






0.0635(04) 


0.0979(09) 


0.75 


0.2183(6) 


0.2509(6) 


0.0668(05) 


0.1028(10) 


0.70 






0.0722(05) 


0.1052(11) 


0.65 


0.2547(6) 


0.2837(6) 


0.0789(06) 


0.1139(10) 


0.60 






0.0868(07) 


0.1167(13) 


0.55 


0.3064(6) 


0.3275(6) 


0.0957(07) 


0.1276(12) 


0.50 






0.1067(09) 


0.1372(13) 


0.45 


0.3834(6) 


0.3897(6) 


0.1171(10) 


0.1473(13) 


0.40 






0.1337(12) 


0.1627(15) 


0.35 


0.5083(6) 


0.4880(6) 


0.1579(14) 


0.1769(18) 


0.30 






0.1877(19) 


0.2034(17) 


0.25 


0.7512(6) 


0.6618(7) 


0.2339(27) 





Table 4. The results for the site-bond percolation thresholds for the 5-dimensional fee and Kagome 
lattices. Error estimates concerning the last digit(s) are indicated between brackets. 



Pb{Ps 



5d- 

P) 



fee 

Ps{Pb 



P) 



5d-Kaj 
Pb{Ps =p) 



;ome 
Ps{pb ^p) 



0.95 


0.0274(5) 


0.0434(5) 


0.1374(16) 


0.2173(20) 


0.85 


0.0319(5) 


0.0474(6) 


0.1555(17) 


0.2284(19) 


0.75 


0.0365(5) 


0.0522(6) 


0.1851(19) 


0.2439(20) 


0.65 


0.0424(5) 


0.0598(6) 


0.2111(19) 


0.2633(19) 


0.55 


0.0500(5) 


0.0661(6) 


0.2538(18) 


0.2957(19) 


0.45 


0.0627(6) 


0.0787(6) 


0.3245(20) 


0.3494(20) 


0.35 


0.0816(6) 


0.0963(6) 


0.4506(20) 


0.4243(19) 


0.25 


0.1179(6) 


0.1283(6) 
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each lattice the percolation threshold was calculated for at least three values of the 
linear lattice size L. The percolation threshold Pc{oo) for infinite lattices can be 
found by fitting these results for different lattice sizes to the scaling relation 

be(L)-pe(00)|(xL-l/^ (4) 

where the critical exponent ly has the value 4/3 in two dimensions, and 0.88, 0.68, 
0.57 in 3, 4, and 5 dimensions respectivelyllJ. These steps are precisely the same 
as for pure site or pure bond percolation. The results are shown in Tables y-^ and 
in Figures ^-^- We have checked that for pure site percolation, and for pure bond 
percolation, our program yields results equivalent to more precise literature values 
(see Table |). 

We verified several scaling laws in percolation, as a check on the program. First, 
we calculated the correlation length ^ below the percolation threshold, as a function 
of pc — p. The correlation length scales as 

^ocipc-py". (5) 

The calculation was performed by first computing the 'gyration radius' R^ for a 
cluster of s sites. This radius is defined by 



R 



2 



1 "" 1 " 

-^|ri-rop, with ro = -^r„ (6) 

i=l i=l 



where r^ is the position of the z*'^ site in the cluster. The correlation length ^ can 
then be calculated as (see Ref.E3, p. 50) 

e = %4^, (7) 

where Ug is the mean number per lattice site of clusters of s sites. We verified 
numerically that the correlation length obeys the scaling relation (|^). Our result 
for site-bond percolation on the Kagome and diamond lattices is shown in Figure g. 
A fit of the scaling relation to the data yields ly = 1.31 ±0.01 for the Kagome lattice, 
and ly = 0.91 ± 0.04 for the diamond lattice. This is within one or two standard 
deviations of the established values oi v — 4/3 in two dimensions and v = 0.88 in 
three dimensions. 

Also we have checked the scaling of the mean cluster size S and the strength of 
the infinite cluster, both as a function of \p — Pc\. The average cluster size S scales 
with the same exponent both above and below the percolation threshold: 

s^\p,-pr, (8) 

with 7 = 43/18 in two dimensions, and 1.80 in threeE2l. The results for site-bond 
percolation on the Kagome and diamond lattices is shown in Fig. |^. The fits to the 
scaling relation (||) yield values for 7 that are consistent with the literature values. 
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diamond lattice: 

fit to ^ oc (p^-p)-^, V = 0.91 ±0.04 
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Fig. 2. The correlation length f as a function of pc — p. The runs were performed at {pt = 
0.65, ps = 0.6284(5)} for the diamond lattice and {pt = 0.75, ps = 0.7757(8)} for the Kagome 
lattice. 
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Fig. 3. The mean cluster size S above and below the percolation threshold, as a function of 
\pc — p\- The runs were performed at {pi, = 0.65, ps = 0.6284(5)} for the diamond lattice (solid 
lines) and {pi, = 0.75, ps = 0.7757(8)} for the Kagomc lattice (dashed lines). 
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The strength P of the infinite cluster scales as 

Pa(p-p,)^ (9) 

with (3 = 5/36 in two dimensions, and 0.41 in threeEJ. The results for site-bond 
percolation on the Kagome and diamond lattices is shown in Fig. H. The fits to the 
scaling relation (pi) yield values for [3 that are consistent with the literature values. 
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Fig. 4. The strength P of the infinite cluster as a function of p - 
^t \_Pb = 0.65, ps = 0.6284(5)} for the diamond lattice and {p(, : 
Kagome lattice. 



Pc- The runs were performed 
: 0.75, Ps = 0.7757(8)} for the 



3. Approximate equations 

It is clear from Figure m that there are systematic deviations from the Yanuka- 
Englman line. The accuracy of the numerical values, as indicated in Tables wM is 
easily high enough to warrant this conclusion. The Yanuka-Englman line is situated 
above the numerical values for all calculated examples. One can say it is an upper 
estimation for the real critical curve. In general it seems that the deviations are 
largest for lattices with a low ratio of a = Pcb/Pcs, i-e. lattices for which the 
bond percolation thresholds is much lower than the site. The most clear example 
is the 'low-ph' lattice (see Fig. ||), which has Pcs — 0.5 and pct — 0.2006(3) (see 
Table H). This lattice has been constructed from the triangular lattice in two steps, 
see Fig. |^. All the lattices in this figure are fully triangulated and have therefore 
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Pes = ^, as shown by Sykes and Essamtj. The bond percolation threshold, however, 
will decrease with every step in Fig. pi and these steps can be repeated even further. 
In this way one can construct lattices with arbitrarily low ratio Pct/Pcs- 




Fig. 5. The triangular lattice (left) can be modified by adding a site in the middle of each triangle 
(middle). The same operation once more yields the lattice on the right. All these lattice have 
Ps = i but the bond percolation threshold decreases from left to right. 

One of the most remarkable results in our calculations is the results for two 
pairs of lattices with very close or equal critical probabilities. One pair is the face 
centered cubic (fee) lattice and the hexagonal close packed (hep) lattice. Another 
pair, in two dimensions, is the direct lattice and the puzzle lattice (we use here the 
nicknames as in Ref. E3). The fee and hep lattices, although they are topologically 
different, have pjl^"^ = p^,^^ and p^^'^ = pjjf^ within the present small error margins 
(see Table ^) . Also the critical curves we calculated for these lattices are the same 
within our error margin. In two dimensions, the pure site and pure bond percolation 
thresholds of the direct and the puzzle lattice are almost the same. Also for this 
pair the critical curves lie close together. 

These results suggest that once the pure site and pure bond thresholds are 
known, the whole critical curve is fixed. In other words, there is no need to use 
extra parameters (such as critical exponents, dimension, or coordination number) 
to describe a critical curve but Pcb and Pcs- In this context it is interesting to note 
that in the limiting case of Bethe lattices, Eq. (g), there are also no parameters 
other than pc- 

One can try to find an equation for the critical curve by making use of the work 
of Hammersley, who derived boundaries for mixed percolationc 

pTiPb) = — , 
Pb 

and 

pf-"(pfc) = ^. 
Pb 

Let us look for a critical relationship as a combination of these boundaries 

pM) = Ci(pb) •^^""(pb) + C2{Pb) -pTiPb)- 

Take the functions C'i{pb) and C2{pb) in the simplest form 

^ / X 1 - Pfc / 1 ior pb^pcb 
^'^""^ ' T^Wb - { iorpb^l ' 
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and 



C2(P5 



Pb - Pcb 
1 -Pcb 

After some transformation we have 



1 for Pb = I 
for Pb = Pcb 



Pb [ps + A) = B, 



where 



A 



B = p. 



Pcb - Pes 
1 - Pcfc 
1 - Pet 



'cb 



1 -Pcb 



(10) 

(11) 

(12) 



The relationship Eq. (|o|) gives almost the same results as Eq. (g) provided ph 
and Ps are close (square, simple cubic, octagonal, etc.). In these cases both curves 
and the numerical points are in good agreement with one another. In other cases 
Eq. (no) gives much better approximations than Eq. (0) (see Figure 0) . 
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Fig. 6. The percolation threshold pairs {ps,Pb) fcr rnany different lattices in various dimensions. 
The points are from the simulations (see Tables lll-W) , the lines are from Eq. ( |lO| ) . 

We examined a number of different relationships with 2 and 3 parameters to 
describe the numerical results. Although the results were excellent for some lattices, 
we could not find a relationship suitable for all lattices. 
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Hence the problem to find the real law for critical probabilities is still unsolved. 
We hope this work with the large collection of high precision numerical results and 
the new estimation of critical probabilities will assist to solve the problem. 
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Table 5. The site and bond percolation thresholds used in the article, together with the references 
to the literature. 



lattice 



Ps [Ref.] 



honeycomb 

Kagome 

square 

triangular 

octagonal 

star-dual 

low-pf, 

puzzle 

direct 

simple cubic 

bcc 

fee 

hep 

diamond 

4d-cubic 
4d-fcc 

5d-fcc 
5d-Kagome 



0.697043(2) 

0.652 703... 

0.592 746 0(5) 

0.5 

0.5 

0.5 

0.5 

0.550806(2) 

0.550213(2) 

0.3116080(3) 
0.245 9615(10) 
0.199 2365(10) 
0.199 24(5) 
0.4301(2) 

0.196901(5) 
0.0842(3) 

0.0431(3) 
0.2084(4) 



Pb [Ref.] 



0.652 703... 

0.524405 3(3) 

0.5 

0.347296... 

0.3230(2) 

0.2593(2) 

0.2006(3) 

0.4142(2) 

0.4195(2) 

0.2488126(5) 
0.180 2875(10) 
0.1201635(10) 
0.12015(5) 
0.3893(2) 

0.1600(2) 
0.049(1) 

0.026(2) 
0.130(2) 



1] 



17,[] 
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